Human Molecular Genetics, 2014, Vol. 23, No. 16 4315-4327 

doi.10. 1093/hmg/ddul48 

Advance Access published on April 4, 2014 

Activating mutations in RRAS underlie a phenotype 
within the RASopathy spectrum and contribute 
to leukaemogenesis 
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RASopathies, a family of disorders characterized by cardiac defects, defective growth, facial dysmorphism, vari- 
able cognitive deficits and predisposition to certain malignancies, are caused by constitutional dysregulation of 
RAS signalling predominantly through the RAF/MEK/ERK (MAPK) cascade. We report on two germline muta- 
tions (p.Gly39dup and p.Val55Met) in RRAS, a gene encoding a small monomeric GTPase controlling cell adhe- 
sion, spreading and migration, underlying a rare (2 subjects among 504 individuals analysed) and variable 
phenotype with features partially overlapping Noonan syndrome, the most common RASopathy. We also iden- 
tified somatic RRAS mutations (p.Gly39dup and p.Gln87Leu) in 2 of 1 1 0 cases of non-syndromic juvenile mye- 
lomonocytic leukaemia, a childhood myeloproliferative/myelodysplastic disease caused by upregulated RAS 
signalling, defining an atypical form of this haematological disorder rapidly progressing to acute myeloid leu- 
kaemia. Two of the three identified mutations affected known oncogenic hotspots of RAS genes and conferred 
variably enhanced RRAS function and stimulus-dependent MAPK activation. Expression of an RRAS mutant 
homolog in Caenorhabditis elegans enhanced RAS signalling and engendered protruding vulva, a phenotype 
previously linked to the RASopathy-causing SHOC2 S2G mutant. Overall, these findings provide evidence of a 
functional link between RRAS and MAPK signalling and reveal an unpredicted role of enhanced RRAS function 
in human disease. 



INTRODUCTION 

Signalling elicited by activated cell surface receptors and trans- 
duced through RAS proteins to the RAF/MEK/ERK and PI3K7 
AKT cascades is central to cell proliferation, survival, differen- 
tiation and metabolism ( 1 ,2). Owing to this nodal role, enhanced 
traffic through RAS proteins and their downstream effectors has 
been established to have a major impact on oncogenesis (3,4). 
This signalling network also controls early and late developmen- 
tal processes (e.g. organogenesis, morphology determination, 
synaptic plasticity and growth), and germline mutations in a 
number of genes encoding transducers and modulatory proteins 
participating in the RAS/MAPK signalling pathway have been 
causally linked to Noonan syndrome (NS) (5), one of the most 
common diseases affecting development and growth, and a 
group of clinically related syndromes, the so-called RASopa- 
thies (6-8). In this family of disorders, constitutional dysregula- 
tion of RAS signalling can be caused by enhanced activation of 
HRAS, KRAS and NRAS (RAS proteins hereafter), aberrant 
function of upstream signal transducers or effectors (PTPN1 1/ 
SHP2, SOS1, SHOC2, RAF1, BRAF, MAP2K1/MEK1 and 
MAP2K2/MEK2) or inefficient down modulation by feedback 
mechanisms (CBL, NF 1 and SPRED 1 ). More recently, RIT1 , en- 
coding a monomeric GTPase structurally linked to RAS pro- 
teins, was identified as disease gene implicated in NS (9), 
extending the concept of 'RASopathy gene' to a transducer 
that contributes to signal propagation through RAS effector 
pathways but does not belong to the RAS/MAPK signalling 
backbone. 

Clinical manifestations of RASopathies include postnatal 
reduced growth, a wide spectrum of cardiac defects, facial dys- 
morphism, ectodermal and skeletal anomalies and variable 



cognitive deficits (5,8,10). Consistent with the key role of most 
RASopathy genes in oncogenesis, these disorders are also char- 
acterized by variably increased risk for certain haematologic ma- 
lignancies and other paediatric cancers (6,7,11,12). Most of 
these conditions are genetically heterogeneous, and the under- 
lying disease gene has not been identified yet for a still significant 
fraction of cases. Based on the strict mechanistic link between 
the molecular events controlling development and contributing 
to oncogenesis, these 'missing' genes represent excellent candi- 
date oncogenes/tumour suppressors. 

Here, we report that constitutional dysregulation of RRAS 
function is associated with a Mendelian trait within the RASopa- 
thy spectrum and that somatically acquired mutations in the 
same gene occur in an aggressive form of juvenile myelomono- 
cytic leukaemia (JMML), a rare childhood myeloproliferative/ 
myelodysplastic neoplasm representing the archetypal somatic 
RASopathy (13), rapidly progressing to acute myeloid leukae- 
mia (AML). We also demonstrate that RASopathy-causing 
RRAS mutations are activating and promote signalling perturb- 
ation by enhancing stimulus-dependent MEK, ERK and, at a 
lower extent, AKT phosphorylation. 

RESULTS 

Identification of candidate disease genes and RRAS mutation 
analysis 

While the core of the machinery implicated in RAS signalling 
has been characterized comprehensively, signal propagation 
through this network is likely to include a larger number of pro- 
teins playing a modulatory or structural role (14), whose aberrant 
or defective function is expected to perturb development and 
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contribute to oncogenesis. Based on this assumption, we used a 
protein interaction/functional association network analysis to 
select a panel of genes encoding proteins functionally linked to 
the RAS signalling network as candidates for NS or a related 
RASopathy (15). Candidate gene selection was based on the 
use of the previously identified RASopathy genes as 'seed' pro- 
teins (i.e. proteins used to build the interaction/functional net- 
works), and considering a panel of databases to construct 
functional subnetworks (Supplementary material, Table S 1 and 
Fig. SI). Sequence scanning of the best candidates in a RASopa- 
thy cohort including 96 unrelated subjects negative for mutations 
in known disease genes allowed the identification of a functionally 
relevant RRAS change (c.l63G>A, p.Val55Met) (Supplemen- 
tary material, Fig. S2) in an adult subj ect with clinical features sug- 
gestive of NS but lacking sufficient characteristics to allow a 
definitive diagnosis (Supplementary material, Table S2). Parental 
DNA was not available for segregation analysis. The mutation 
was not identified among >400 population-matched unaffected 
individuals, indicating that it did not represent a common poly- 
morphic nucleotide substitution. This change, rs368625677 
(dbSNP 138), had been described in 1/13,006 alleles in the 
NHLBI Exome Sequencing Project (http://eversusgs.wa 
shington.edu/EVS/). Of note, similar frequencies have been 
reported in the same database for recurrent RASopathy-causing 
mutations (e.g. c.922A>G in PTPN11, and c.1259G>A in 
CBL). Mutation analysis was extended to additional 408 patients 
with NS or a clinically related phenotype tested negative for muta- 
tions in the major NS disease genes (see Materials and Methods), 
allowing to identify one sporadic case heterozygous for a three- 
nucleotide duplication (c.H6_118dup, p.Gly39dup) (Supple- 
mentary material, Fig. S2). Parental DNA sequencing of the 
relevant exon demonstrated the de novo origin of the variant, 
and STR genotyping confirmed paternity. In this subject, the du- 
plication was documented in DNA obtained from skin fibroblasts, 
excluding a somatic event restricted to haematopoietic cells. The 
subject had features reminiscent of NS (Fig. 1 A and Supplemen- 
tary material, Table S2), with onset of AML suspected to represent 
a blast crisis of JMML (Supplementary material, Table S3 and 
Fig. S3). In this patient, exome sequencing performed on leu- 
kaemic and non-leukaemic DNA failed to disclose any additional 
relevant germline/somatic change affecting genes known to be 
mutated in RASopathies and JMML, as well as genes directly 
linked to the RAS signalling network, further supporting the 
causal role of the identified RRAS lesion. Based on this associ- 
ation, the occurrence of RRAS mutations was also explored in a 
panel of genomic DNAs obtained from bone marrow aspirates/ 
circulating leukocytes of 110 subjects with JMML. Heterozy- 
gosity for the previously identified Gly 39 duplication and the 
c.260A>T (p.Gln87Leu) change was observed in two patients 
with JMML rapidly progressing to AML (Supplementary mater- 
ial, Table S3 and Fig. S3). Both lesions were absent in non- 
leukaemic DNA, indicating their somatic origin (Supplementary 
material, Fig. S2). These subjects also carried a somatic NRASmu- 
tation, suggesting that the two hits might cooperate with this par- 
ticularly severe form of disease. Sequencing of isolated JMML 
myeloid colonies in patient 14385 showed that NRAS and RRAS 
mutations coexisted in the same progenitors but failed to establish 
their sequence of appearance during leukaemogenesis, not allow- 
ing to discriminate whether the latter was involved in initiation or 
progression of disease. 



Structural analyses 

RRAS encodes a 23-kD a membrane -bound monomeric GTPase 
with 55-60% amino acid identity to RAS proteins (16). This 
highly conserved structure is flanked by a unique 26-amino 
acid region at the A^-terminus (Fig. IB). Similarly to the other 
RAS family proteins, RRAS binds to GTP and GDP with high 
affinity and specificity and functions as a molecular switch by 
cycling between active, GTP-bound and inactive, GDP-bound 
states (17). RRAS is activated by guanine nucleotide exchange 
factors (GEFs) in response to signals elicited by cell surface 
receptors. In the GTP-bound state, two functionally conserved 
regions, switch I and switch II (Fig. IB), undergo a conform- 
ational change enabling RRAS to bind to and activate effector 
proteins. This interaction is terminated by hydrolysis of GTP 
to GDP, which is promoted by GTPase-activating proteins 
(GAPs) and results in switching towards the inactive conform- 
ation. Disease-associated RRAS mutations affected residues 
highly conserved among orthologs andparalogs (Supplementary 
material, Fig. S4) residing in the GTP-binding pocket (Fig. 1C) 
and were predicted to be damaging with high confidence (Sup- 
plementary material, Table S4). Among them, Gin 87 , homolog 
of Gin 61 in RAS proteins, is directly involved in catalysis 
(18,19). The p.Gln87Leu substitution had previously been 
reported as a rare somatic event in lung carcinoma, and mutations 
affecting Gin 61 are among the most recurrent oncogenic lesions in 
RAS genes (COSMIC database, http://cancer.sanger.ac.ulc/ca 
ncergenome/projects/cosmic/). Likewise, p.Gly39dup altered 
the Gl motif participating in GTP/GDP binding and GTPase ac- 
tivity (Fig. IB). Within this motif, Gly 12 and Gly 13 (Gly 38 and 
Gly 9 in RRAS) represent major mutation hot-spots in human 
cancer (COSMIC database) and account for the majority of germ- 
line HRAS mutations causing Costello syndrome (20). Of note, 
analogous insertions in RAS proteins have been reported in 
JMML and other malignancies (2 1 -24). In contrast, no somatic/ 
germline RAS mutation affecting Val 29 , homolog of Val 55 in 
RRAS, had previously been reported. Val 55 side-chain is not dir- 
ectly involved in GTP/GDP binding, GTP hydrolysis or inter- 
action with effectors. However, it has been reported that 
H-bonds are possible between the backbone of Val 9 in HRAS 
and GDP/GTP (25). Furthermore, it has been suggested that 
Val 29 can play a role in the transition between the GDP- and 
GTP-bound states (26), as supported by the evidence that the 
Val29Gly substitution in HRAS accelerates the GDP/GTP ex- 
change in vitro (27). 

Molecular dynamics (MD) simulations were performed to 
predict in silico the effects of p.Val55Met on the structure and 
dynamics of RRAS (Fig. 2). The mutation was introduced in 
the available crystallographic structure of RRAS in complex 
with GDP and Mg 2+ , and the system was simulated in water 
for 200 ns. For comparison, MD simulations were also per- 
formed using the wild-type protein, which maintained a stable 
structure along the whole simulation, as expected (Fig. 2 A, left 
panel). In contrast, a dramatic local structural transition extend- 
ing up to the switch I region (residues 58-64), which mediates 
effector binding, was documented for the RRAS V55M mutant, 
after ~80 ns (Fig. 2A, right panel). This conformational transi- 
tion resulted in an increased solvent exposure of Met 55 , in agree- 
ment with the higher hydrophilicity of this residue compared 
with Val, and was accompanied by the formation of a stable 
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Figure 1. RASopathy-causing and leukaemia-associated RRAS mutations. (A) Facial features of the affected subject (9802) heterozygous for the de novo germline 
c. 1 1 6_1 1 8dup. (B) RRAS exon- intron arrangement with coding exons as blue boxes. RRAS functional motifs include the GTP/GDP binding domain (G 1 to G5 , start- 
ing from the iV-terminus) (red), switch I (light green), switch II (dark green) and hypervariable region (light brown) with the C-terminal CAAX motif (dark brown). The 
unique JV- terminal region is also shown (violet). Location of disease-associated mutations is reported. (C) Position of affected residues on the three-dimensional struc- 
ture of RRAS in its GDP-bound, inactive state (PDB: 2FN4) (above) and that of non-hydrolysable GTP analogue (GppNHp)-bound, active HRAS (PDB: 5P21) 
(below). The red surface indicates Gly 39 and Val 40 (Gly 13 and Val 14 , in HRAS), whereas Val 55 (Val 29 ) and Gin 87 (Gin 61 ) are shown in blue and green, respectively. 
GDP is reported as semi-transparent yellow surface. 



cluster involving residues He 50 , Met 55 and Tyr 58 (Fig. 2A and 
Supplementary material, Table S5) permitted by the unbranched 
and long side-chain of Met 55 . No further significant conform- 
ational changes were observed for the remaining interval of 
the simulation. The major effect of this structural rearrangement 
was to increase exposure of GDP to the solvent (Fig. 2B), with an 
almost doubled solvent accessible surface area of the nucleotide 
after the conformational transition. This structural rearrange- 
ment was accompanied by a perturbation of the intermolecular 
H-bond network stabilizing GDP binding, with loss of the 
H-bonds between residues at codons 55 and 56, and GDP (Sup- 
plementary material, Table S5). Of note, apossible impact of the 
described structural transition on RRAS binding to GEF pro- 
teins, which bind to this region and mediate GDP release, was 
also noticed. Specifically, we observed that after the conform- 
ational rearrangement, the RPvAS 5 region implicated in 
GEF binding populated a structure similar to that assumed in 
RAS/GEF complexes (Fig. 2C), suggesting a possible enhanced 
interaction of the disease-associated RRAS mutant with GEFs. 
In particular, Tyr 58 was observed to adopt a side-chain orienta- 
tion very similar to that of the RAS homolog Tyr 32 in the 
HRAS/SOS1 complex, which has been shown to contribute to 
the structural rearrangements of switch I and interaction with 
GEFs (28-30). 

Overall, these data supported an activating role of 
p.Val55Met through enhanced GDP release as a result of a 
decreased affinity for the nucleotide and/or enhanced inter- 
action with a GEF. 



Biochemical and functional characterization 
of RRAS mutants 

Previous studies documented the gain-of-function role of 
p.Gln87Leu on RRAS function, and MAPK and PBK/AKT sig- 
nalling (31). To characterize the impact of p.Val55Met and 
p.Gly39dup on protein function, we analysed the intrinsic and 
GEF-accelerated nucleotide exchange reaction of these mutants. 
Dissociation kinetics analysis demonstrated a dramatically in- 
creased intrinsic (RRAS G35idup ) and GEF-stimulated (RRAS G39dup 
and RRAS V55M ) dissociation rate of mantGDP, indicating a facili- 
tated nucleotide release in both mutants (Fig. 3A). RAS proteins 
exhibit low intrinsic GTPase activity that is enhanced by GAPs 
(32). Assessment of RRAS G39dup and RRAS V55M GTPase activity 
documented a significantly reduced intrinsic and GAP-stimulated 
GTP hydrolysis in the former (Fig. 3B and Supplementary mater- 
ial, Fig. S5). Finally, the interaction of RRAS proteins with various 
effectors was analysed by fluorescence polarization (Fig. 3C). 
While RRAS WT was found to bind to RAF1, RALGDS, 
RASSF5 and PLCE1 less efficiently than HRAS, an increased 
binding affinity to PIK3CA was observed. Compared with 
RRAS T , aberrant binding behaviour of the two RRAS mutants 
was demonstrated, with RRAS G39dup exhibiting an increased 
binding affinity towards PIK3CA, RAF1, PLCE1 and RASSF5, 
and RRAS V55M to RALGDS. 

To gain further insights into the impact of disease-causing 
mutations on RRAS functional dysregulation and explore their 
effects on RAS signalling, the activation state of RRAS proteins 
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Figure 2. Molecular dynamics (MD) simulations. (A) Structural perturbations promoted by the p.Val55Met substitution as obtained from MD simulations of the 
RRAS/GDP complex. The wild-type (WT) protein is also shown for comparison. Top panels report the protein structures at the beginning of simulations, whereas 
the final structures (200 ns) are shown at the bottom. The final structure of RRAS V55M is well representative of the last 120 ns of the trajectory. The protein 
surface of RRAS is shown with GDP (yellow). The mutated residues and those forming a cluster in the simulation of mutated RRAS are coloured as follows: 
Val 55 /Met 55 (blue), Tyr 58 (pink) and He 50 (cyan). Residues 59-64, which, together with Tyr 58 , form the switch I region, are coloured in green. (B) Solvent accessible 
surface of GDP in the MD simulations of wild-type (red) and mutant (blue) RRAS/GDP complexes. (C) Conformation of the loop comprised between Val 55 /Met 55 and 
Asp 59 in wild-type (red) and mutant (blue) RRAS/GDP complexes obtained from MD simulations. GDP is reported as semi-transparent yellow surface. Superimposed 
conformations of the corresponding loop (residues 29-33) in GDP-bound HRAS (violet) (PDB: 4Q21) and GDP-bound HRAS complexed with SOS1 (cyan) (PDB: 
1BKD) are shown for comparison. The side chains of Tyr 58 and the corresponding residue in HRAS, Tyr 32 , are displayed as sticks. 



and extent of signalling through the MAPK and PI3K/AKT 
cascades were evaluated using transient expression in COS-7 
cells. Consistent with the above-mentioned findings, pull-down 
assays revealed a variably higher proportion of active, GTP- 
bound form for both mutants (Fig. 4A). Moreover, similarly to 
what observed under cell-free conditions, RRAS G39dup was re- 
sistant to GAP stimulation. Expression of both mutants promoted 
enhanced serum-dependent MEK, ERK and AKT phosphoryl- 
ation (Fig. 4B), which was more evident in cells expressing the 
RRAS G3?dup mutant. 

Caenorhabditis elegans studies 

To explore further the functional impact of the RASopathy 
causative RRAS mutants on RAS signalling in vivo, we used 
the nematode C. elegans as an experimental model. In C. 
elegans, the role of ras-1, the RRAS ortholog (33), has not 
been characterized yet. On the contrary, proper signalling 
through LET-60, the C. elegans ortholog of RAS proteins, has 
been established to play a crucial role in vulval development 
(34). In particular, LET-60/RAS is known to mediate the 
priming signal (LIN-3/EGF) released by the anchor cell to 
induce the three nearby vulval precursor cells (VPCs), P5.p, 
P6.p and P7.p, to generate a normal vulva. Enhanced and 
decreased signalling through LET-60 and the MAPK cassette 



results in multiple ectopic pseudo vulvae (multi vulva phenotype) 
and a failure in VPC induction (vulvaless phenotype), respect- 
ively (34,35). 

Multiple transgenic lines were generated to conditionally 
express the wild-type ras-1 cDNA (ras-l WT ) or the allele hom- 
ologous to the disease-associated three-nucleotide duplication 
(ras-l G27dup ), which was identified to occur both as a germline 
and somatic event. Exogenous RAS-1 expression was induced 
by heat shock at early L3 larval stage to investigate the effects 
of the mutant protein on vulval development. Animals expres- 
sing ras-1 2 up displayed abnormal vulval morphogenesis 
resulting in the formation of a protruding vulva (Pvl) (Fig. 5A 
and B and Supplementary Material, Table S6), a phenotype asso- 
ciated with aberrant traffic through different signalling cascades 
(36,37). Of note, this phenotype had previously been reported in 
worms expressing the RASopathy causative SHOC2 S2 mutant 
(38). Like those animals, ras-l G 7dup worms showed decreased 
egg-laying efficiency (Egl phenotype), and accumulation of 
larvae inside the mother (Bag-of-worms phenotype). A signifi- 
cantly less penetrant phenotype was observed in animals expres- 
sing ras-1 T . These findings, together with the observation that 
animals lacking ras-1 do not exhibit any vulval defect (Worm- 
Base, http://www.wormbase.org/, and our personal assessment), 
supported the gain-of-function role of the mutation on 
RAS-1 function. At the late L3/early L4 larval stage, vulva 
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Figure 3. In vitro biochemical characterization of the RRAS G39dup and RRAS V55M mutants. (A) Intrinsic (left) and SOS 1 -accelerated (right) mantGDP nucleotide 
dissociation measured in the presence of 20-fold excess of non-labelled GDP. The decrease in mant fluorescence was fitted by single exponentials to obtain k obs values 
for nucleotide dissociation. RRAS G39dup exhibited a 35-fold increased intrinsic dissociation of mantGDP, whereas SOS 1 -accelerated mantGDP dissociation was aug- 
mented by ~ 30-fold for both mutants, compared with RRAS WT and HRAS WT . The values are an average of five to seven independent measurements. (B) Intrinsic 
GTP hydrolysis kinetics (left) and rate constants (right) of RRAS G39dllp and RRAS V55M proteins, documenting the impaired catalytic activity in the former. The ¥^, bs 
values are an average of five to seven independent measurements. (C) Binding of RRAS T , RRAS V55M and RRAS G3 dup to the RAS-binding domain of RAF 1 mea- 
sured as variation in fluorescence polarization of each mantGppNHp-bound RRAS protein at increasing concentrations of RAF1-RBD (left), and dissociation con- 
stants (Kj) for the interaction of HRAS WT and RRAS proteins to the RBDs of RAF 1 , RALGDS, PLCE1 , PIK3CA and RASSF5 (right). Kj values were obtained by 
quadratic fitting of the concentration-dependent binding curves from fluorescence polarization measurement as exemplified for RAF 1 -RBD. Of note, RRAS WT binds 
to RAF1 , RALGDS, RASSF5 and PLCE1 less efficiently than HRAS, whereas an increased binding affinity to PIK3CA is observed. 



morphogenesis normally begins with the descendants of VPC 
P6.p detaching from the cuticle and forming a symmetric invagin- 
ation (Fig. 5C) (34). Animals in which the expression of ras-1 
had been induced at early L3 largely maintained this pattern (17/ 
20). In contrast, in larvae expressing ras-l G27d " p , descendants of 
VPCs P5.p and/or P7.p more frequently detached from the 
cuticle, resulting in larger and more asymmetric invaginations 
(10/30). This morphogenesis defect was the earliest detectable 
effect of the ras-1 7 up allele on vulval development, similarly 
to that previously documented in transgenic lines expressing 
SHOC2 fe2G (38). 

Genetic interaction between the RAS-1 /RRAS mutant and 
LET-60/RAS was also investigated. While expression of the 
RAS-1 G27dup mutant was able to exacerbate the multivulva 



phenotype associated with a hyperactive let-60 allele (nl046), ex- 
pression of wild-type RAS-1 failed to do so (Table 1). Similarly, a 
significant, although partial rescue of the VPC induction defect asso- 
ciated with a let-23/EGFR hypomorphic allele (syl) was observed 
in animals expressing the activating RAS-1 G27dup mutant, but not 
in worms expressing the wild-type counterpart (Table 1). Overall, 
these experiments provided evidence of a positive modulatory 
role of the RAS-1/RRAS mutant on LET-60/RAS signalling. 

DISCUSSION 

Mutations of genes coding for proteins with role in RAS signal- 
ling and the RAF/MEK/ERK cascade have been identified as the 
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Figure 4. RRAS G39dup and RRAS V55M signalling activities in cells. (A) Determination of GTP-bound RRAS levels in COS-7 cells transiently expressing wild-type or 
mutant FLAG-tagged RRAS proteins. Assays were performed in the presence of serum (above), and in serum-free conditions ( — GAP) or in the presence of the neu- 
rofibromin GAP domain (+ GAP) (below). RRAS G39dup was predominantly present in the active GTP-bound form and was resistant to GAP stimulation, whereas a 
slightly increased level of GTP-bound RRAS V55M was observed in the presence of serum. Representative blots of three performed experiments are shown. (B) De- 
termination of MEK, ERK and AKT phosphorylation levels (pMEK, pERK and pAKT) in transiently transfected COS-7 cells cultured in medium with serum (left) or 
basal medium (right). Expression of each RRAS mutant resulted in variably enhanced MEK, ERK and also partially AKT phosphorylation after stimulation. Total 
MEK, ERK and AKT in cell lysates are shown for equal protein expression and loading. Expression levels of exogenous, FLAG-tagged RRAS in cell lysates are shown 
for each experiment. Representative blots of three performed experiments are shown. 



molecular cause underlying a group of clinically related devel- 
opmental disorders, the RASopathies. Here, we used a gene can- 
didacy approach based on large-scale protein-protein 
interaction/functional network analysis to identify RRAS as a 
novel gene implicated in a condition with features within the 
RASopathy spectrum. Disease-causing RRAS mutations are ac- 
tivating and act by maintaining the GTPase in its GTP-bound 
active state. Aberrant RRAS function was demonstrated to 
perturb variably intracellular signal flow through the RAF/ 
MEK/ERK cascade, and to a certain extent also the PI3K7AKT 
pathway. Of note, these gain-of-function mutations are likely 
to define a novel leukaemia-prone condition. Consistent with 
this view, the same class of RRAS lesions was identified to 
occur as acquired somatic event in JMML, characterizing a 
subset of this myeloproliferative/myelodysplastic disorder 
with rapid progression to AML. 

RRAS shares several biochemical properties with HRAS, 
NRAS and KRAS, as well as some common function, including 



stimulation of cell proliferation, survival and transformation 
(19,39). Despite these similarities, however, previous observa- 
tions have emphasized the role of RRAS in cell adhesion, spread- 
ing and migration, and its modulatory function on effectors 
distinct from those used by 'classical' RAS proteins (40,41). 
While PI3R7AKT has been recognized as a major effector 
pathway of RRAS, only a minor impact on MAPK signalling 
had been reported (41,42). The present in vitro findings 
provide evidence that disease-associated RRAS mutants 
enhance the activation of the MAPK cascade, at least in response 
to specific stimuli. On the other hand, the identification of RRAS 
as a novel disease gene implicated in a RASopathy disorder 
further emphasizes the relevance of dysregulated signalling con- 
trolling cell spreading and migration in certain features of NS 
(e.g. congenital heart defects and lymphedema) and JMML 
(leukocyte infiltration in non-haematopoietic tissues) (43-45). 

Caenorhabditis elegans studies provided evidence for a 
genetic interaction between the RAS-l G27dup /RRAS G39dup and 
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Figure 5. Consequences of ras-l a2? up expression on C. elegans vulval development. (A) Heat-shock-driven expression of ras-l WT &ndras-l c27d " p at early L3 stage 
results in protruding vulva (Pvl), egg laying defective (Egl) and bag-of-worms (Bag) phenotypes. Isogenic animals that had lost the transgene (control group) and 
worms expressing the heat shock-inducible vector (empty vector) were subjected to heat shock and scored in parallel for comparison. The dose at which the transgene 
has been inj ected is reported at the bottom. Error bars indicate SD of three independent experiments. Asterisks indicate significant differences compared with ras-1 m 
at the corresponding dose of injection (*P < 0.05; **P < 0.005; ***p < 0.0005; Fisher's Exact Test). (B) A proper vulva develops in heat-shocked control animals 
(left), whereas a protruding vulva is observed in heat-shocked ras- l a27d " p young adults (middle) and adult worms (right). (C) Nomarski images of vulval precursor 
cells in late L3 (left), early L4 (middle) and mid-late L4 (right) stages from synchronized animals heat-shocked at early L3. In control animals (N = 48), only P6.p 
descendants invaginate (upper panel), whereas in 10 of 30 analysed ras-l a2? u/> -expressing worms, P5.p and/or P7.p descendants also detach from the cuticle, gen- 
erating asymmetric invaginations (lower panel) . Black arrowheads point to P6.p descendant invagination, whereas white arrowheads point to P5 .p and P7.p descendant 
invagination. Anterior is to the left and dorsal is up, in all images. 



LET-60/RAS in vivo. Specifically, expression of the RAS-1 
mutant protein was able to rescue, in part, the VPC induction 
defect resulting from a hypomorphic LET-23 mutant and 
enhanced the multivulva phenotype associated with a LET-60 
gain-of-function genetic background. No impact of wild-type 
RAS-1/RRAS expression was observed in both models. We 
also observed that worms expressing ras-l G27dup displayed ab- 
normal vulval morphogenesis (protruding vulva), possibly 
resulting from aberrant morphogenetic movements of the VPC 
descendant cells. Of note, we observed an equivalent phenotype 
in transgenic lines expressing SHOC2 S2G (38) and a PTPN1 1/ 
SHP2 gain-of-function mutant (our unpublished data), suggest- 
ing functional equivalence of these mutants. Genetic studies 
support the view that these vulva defects arise, in part, through 
perturbation of signalling mediated by the RHO-related 
GTPase, RAC, which plays a critical role in vulval morphogen- 
esis (37). This finding is inline with the established role of RRAS 
on RAC signalling (40,41) and with preliminary data indicating 
enhanced migration and chemotactic capabilities in cells stably 
expressing the disease-associated RRAS mutants (our unpub- 
lished data). 

The biochemical characterization of disease-associated-ft^S 
mutations provided strong evidence for the existence of distinct 
structural and mechanistic effects resulting in an overall en- 
hancement of RRAS signalling. Function of RAS family pro- 
teins in signal transduction is controlled by two events, the 
GDP/GTP exchange and GTP hydrolysis. Any perturbation of 
these processes can affect dramatically the fine-tuned balance 
of the GTPase interaction with effectors and signal output. The 
majority of gain-of-function mutations affecting RAS proteins, 
including those contributing to oncogenesis, trigger the accumu- 
lation of these GTPases in the active state by impairing intrinsic 
GTPase activity, and/or conferring resistance to GAPs (17). This 
is also the case of two of the three mutations identified in this 



study, p.Gly39dup and p.Gln87Leu, the latter corresponding to 
the p.Gln61Leu in RAS proteins (present study and ref. 18,19). 
The characterization of the biochemical behaviour of 
RRAS G39dup , however, also demonstrated a dramatic increase 
in both the intrinsic and GEF-catalysed nucleotide exchange as 
a process contributing to the accumulation of this mutant in its 
GTP-bound state. Aberrant GEF-accelerated nucleotide ex- 
change dynamics was identified as the event driving functional 
dysregulation in the RRAS V55M mutant, which was documented 
to be hyper responsive to GEF stimulation, but retained stimulus- 
dependency. Remarkably, the RRAS G39dup and RRAS V55M 
mutants were demonstrated to exhibit a diverse binding behav- 
iour to effectors suggesting a differential impact of mutations 
on downstream signalling cascades, including PI3K/AKT and 
RALGDS/RAL, whose biological significance and impact, 
however, require further studies. 

The clinical phenotype of the two subjects with germline 
RRAS mutations was reminiscent of NS. The individual with 
the Gly 39 duplication displayed pulmonic stenosis, reduced 
growth, cafe-au-lait spots, mild motor delay and low-set ears, 
which recur in NS (5). Facial features, however, were distinctive, 
and not typical of NS. In contrast, the patient heterozygous for 
the p.Val55Met substitution exhibited a very mitigated pheno- 
type characterized by suggestive facial characteristics (triangu- 
lar face, downslanting palpebral fissures and low-set ears), low 
posterior hairline, broad chest and borderline cognitive abilities, 
without cardiac involvement or defective growth, indicating that 
clinical features associated with RRAS mutations might be quite 
subtle. Of note, the milder phenotype associated with the 
p.Val55Met change is consistent with the weaker perturbing 
effect of the RRAS V55M mutant on MAPK and PI3K/AKT sig- 
nalling compared with the RRAS G39dup protein. Additional in- 
formation on the spectrum of germline RRAS mutations, their 
associated phenotype and their functional impact on signalling, 
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Table 1. Vulva phenotypes in C. elegans mutant strains expressing wild-type RAS-1 or the disease-associated RAS-1 ' up mutant 



Genotype 


Transgene 


N 


Muv (%) 


Vul (%) 


Pvl (%) 


N 


P6.p induction (%) 


wild-type 


none 


207 


0 


0 


1.0 


48 


100 


let-60{nl()46) 


none 


201 


77.9 




0.5 


50 


100 


let-60(nl()46) 


ras-l™ 


244 


76.4 




2.8 


43 


too 


let-60(nl()46) 


ras-l G27dup 


231 


87.1" 




3.0 


50 


100 


let-23(syl) 


none 


194 




87.8 


3.6 


178 


13.4 


let-23(syl) 


ras-l^ 


169 




84.3 


4.1 


156 


14.0 


let-23(syl) 


ras-l G27dup 


282 




83.3 


10.3° 


128 


24.2° 



Strains: let-60(nl046) is a gain-of-function allele of let-60 (ortholog of the human HRAS, KRAS and NRAS genes); let-23(syl) is a hypomorphic allele of let-23 
(ortholog of the human EGFR gene). ras-l WT an&ras-l G27d " p indicate hsp-16.41::ras-l WT - and hsp-16.41::ras-l '^''"''-containing constructs injected at 100 ng/(jd, 
respectively. After each cross, isogenic worms that had lost the transgene were cloned separately and used as controls. 

Animals were grown at 20 "C and heat-shocked in parallel at early L3 stage. ./Vindicates the number of animals scored. Multivulva (Muv), vulvaless (Vul) and 
protruding vulva (Pvl) phenotypes are expressed as percentage of worms with ectopic pseudovulvae, animals lacking a vulva and adults with a protruding vulva, 
respectively. Induction of vulval cell fate is expressed as percentage of P6.p that has been induced to invaginate. 
In all comparisons, /"-values were calculated using two-tailed Fisher's exact test. 
Significantly different from let-60(nl046) (P < 0.02). 

Significantly different from let-23(syl) (P < 0.01) and let-23 (sy 1) ;ras- 1 WT (P < 0.02). 
Significantly different from let-23(syl) (P = 0.02) and let-23(syl);ras-l WT {P < 0.05). 



however, is necessary to establish clinically relevant genotype - 
phenotype correlations. 

JMML is a clonal myeloproliferative/myelodysplastic dis- 
order of childhood characterized by overproduction of immature 
myeloid cells that variably retain the capacity to differentiate. 
Upregulation of RAS/MAPK signalling owing to germline and 
somatic mutations in PTPN11, NRAS, KRAS, NF1 and CBL is 
a major event implicated in this malignancy (13,46,47). Our 
data document that upregulated RRAS function represents a 
novel event contributing to JMML pathogenesis and/or disease 
progression. Notably, somatic RRAS mutations co-occurred 
with acquired NRAS lesions in atypical JMML characterized 
by late onset and rapid progression to AML. While JMML is gen- 
erally an aggressive malignancy, a subset of NRASIKRAS 
mutation-positive patients has been reported to exhibit a mild 
course, with spontaneous remission despite the T^S-mutated 
clone persisting for years (48-50, our unpublished data). This 
suggests that in some instances, certain NRAS mutations are 
not sufficient to support full leukaemogenesis, requiring syner- 
gism with a second RAS signalling targeting event. In line 
with this view, NRAS mutations have been documented to 
co-exist with defects in other RASopathy genes (e.g. PTPN11) 
in some cases resulting in a particularly aggressive disease re- 
sembling AML with myelodysplasia-related changes (51,52), 
as that observed in the present cases. Other studies, however, 
are required to appreciate more precisely the role of enhanced 
RRAS function in leukaemogenesis as well as its clinical rele- 
vance in haematological malignancies. 

In conclusion, our findings document that germline activating 
mutations in RRAS underlie a condition within the RASopathy 
family that may resemble NS phenotypically. In the examined 
cohorts, RRAS lesions were found to account for only a small 
portion of cases, which might be related to their severe conse- 
quences on embryonic/foetal development and/or to the biased 
selection of the subjects included in this study. Based on the 
present findings, however, RRAS mutations are expected to be 
more common among subjects with clinical features only par- 
tially overlapping NS, and particularly in patients with syndrom- 
ic JMML/ AML not associated with mutations in the PTPN11, 
NF1, CBL, KRAS and NRAS genes. While further efforts are 



required to characterize more precisely the clinical impact of 
germline and somatic mutations affecting RRAS, our findings 
suggest an unpredicted role of this GTPase in development 
and haematopoiesis. Consistent with the recent identification 
of RIT1 as disease gene implicated in a significant proportion 
of NS (9), our findings further extend the concept of 'RAS- 
opathy gene' to a transducer whose dysregulated function per- 
turbs signal flow through the MAPK cascade but does not 
belong to the core RAS/MAPK signalling cassette. 

MATERIALS AND METHODS 

Selection of RASopathy candidate genes 

A web-based tool, Genes2FANs (http://actin.pharm.mssm.edu/ 
genes2FANs), using a large-scale protein -protein interaction 
network coupled to a panel of functional association networks 
(FANs) was utilized to build a subnetwork connecting proteins 
to the known RASopathy genes (i.e. PTPN11, SOS1, NF1, 
SPRED1, CBL, NRAS, KRAS, HRAS, RAF1, BRAF, SHOC2, 
MAP2K1 and MAP2K2), as seed proteins. Gene Ontology (bio- 
logical process tree), mammalian phenotype browser, and Con- 
nectivity Map (drug-associated gene expression signatures), 
ChEA and TRANSFAC (transcription factor networks) data- 
bases were selected to construct the functional subnetworks uti- 
lized for prioritization of candidates (15). The programme 
allows to calculate z-scores for intermediate nodes (i.e. candi- 
dates), which are ranked based on their connections to the seed 
proteins. 

Subjects and mutation analysis 

Three cohorts of patients were considered in the study. A first 
group including 96 subjects with clinical features within the 
RASopathy spectrum and without mutation in previously identi- 
fied RASopathy genes (i.e. CBL, PTPN11, SOS1, KRAS, HRAS, 
NRAS, RAF1, BRAF, SHOC2, MAP2K1 and MAP2K2) was 
screened for a selected panel of candidates. A second cohort in- 
cluding 408 subjects with NS or a closely related phenotype pre- 
viously tested negative for mutations in a heterogeneous subset 
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of RASopathy genes was seamed for RRAS mutations. All sub- 
jects included in this cohort had been screened for mutations in 
PTPN11, SOS1 and RAF1 genes. In both cohorts, the clinical 
diagnosis was made on the basis of standardized clinical criteria 
assessed by experienced clinical geneticists and paediatricians. 
Clinical features for the majority of subjects satisfied diagnostic 
criteria reported for NS, LEOPARD syndrome and cardiofacio- 
cutaneous syndrome (53-57), but individuals who lacked suffi- 
cient features for a definitive diagnosis were also included. RRAS 
mutation analysis was also carried out on a cohort including 110 
subj ects with non-syndromic JMML that had prospectively been 
collected and genotyped (58). Mutation screening was perfor- 
med on the entire RRAS coding sequence and flanking intronic 
stretches (NC_000019.10, 49635295.. 49640143, complement; 
NM_006270.3; NP_006261.1) on genomic DNA extracted 
from circulating leukocytes (cohorts I, II and III) or bone marrow 
aspirates (cohort III) by denaturing high-performance liquid 
chromatography (DHPLC) (3 1 00 or 3 500HT WAVE DNA frag- 
ment analysis system, Transgenomic, Omaha, NE, USA) and/or 
direct bidirectional sequencing (ABI Prism 3 1 30, 3730 and 3500 
Genetic Analyzers; Applied Biosystems, Foster City, CA, USA). 
Primer pairs, PCR and DHPLC conditions are available upon 
request. dbSNP137 (http://www.ncbi.nlm.nih.gov/projects/ 
SNP/snp_summary.cgi), HapMap (rel.27) (http://hapmap.ncbi. 
nlm.nih.gov/) and 1000 Genomes (http://www.l000genomes. 
org/) databases were used to annotate the identified sequence 
variants. SIFT (http://sift.jcvi.org/), PolyPhen-2 (http://genetics. 
bwh.harvard.edu/pph2/) and MutationTaster (http://www.muta 
tiontaster.org/) were used to predict the functional impact of 
the identified variants. Paternity was confirmed by STR geno- 
typing, using the PowerPlex 16 System (Promega). DNA from 
leukocytes, hair bulb cells, bone marrow aspirates and skin fibro- 
blasts was extracted using standard protocols. DNA specimens 
were collected under Institutional Review Board-approved pro- 
tocols. Informed consent for DNA storage and genetic analyses 
was obtained from all subjects. Permission was obtained to pub- 
lish picture of patient 9802, whereas subject NS1166 denied 
consent for picture publication. 

Exome sequencing 

Targeted enrichment and massively parallel sequencing were 
performed on genomic DNA extracted from circulating leuko- 
cytes and fibroblasts of patient 9802. Exome capture was 
carried out using the SureSelect Human All Exon V4+UTRs 
(Agilent), and sequencing with a HiSeq2000 instrument (Illu- 
mina). Image analysis and base calling were performed using 
the Real Time Analysis (RTA) pipeline v. 1.14 (Illumina). 
Paired-end reads alignment to the reference human genome 
(UCSC GRCh37/hgl9) and variant calling were carried out 
using the CASAVA v. 1.8 pipeline (Illumina). Variant annota- 
tion, SNP filtering (dbSNP135, 1000 Genomes, HapMap and 
IntegraGen Exome databases) and patient-matched germline 
variant filtering were attained using an in-house pipeline by Inte- 
graGen (Evry, France). 

Molecular dynamics analyses 

Starting coordinates for MD simulations were obtained from the 
RRAS crystallographic structure in complex with GDP and 



Mg 2+ (PDB: 2FN4; RCSB Protein Data Bank, http://www. 
rcsb.org/pdb/home/home. do). The /V-terminus and C-terminus 
of RRAS, absent in the crystal, were not considered in simula- 
tions. Deep View software was used to add atoms missing in 
the PDB file, belonging to residues 31, 96, 1 14, 121 and 192 as 
well as to substitute residue 55 (simulation of the p.Val55Met 
mutant) (59). All MD simulations were performed with 
GROMACS 4.5 package, by using the GROMOS96 43al 
force field parameters for the protein. Parameters for GDP 
were taken from the GROMACS website (http://www.groma 
cs.org). Simulations were performed as previously described 
(60,61), except for some details. Briefly, the proteins were ini- 
tially placed in a dodecahedral box, solvated with ~6700 SPC 
water molecules, and six Na + ions were added to neutralize 
the protein charge. Following initial energy minimization and 
a 100-ps MD run, during which the protein atoms were position 
restrained, the temperature of solute and solvent was raised from 
50 to 300 K in a stepwise manner. 200-ns-long simulations were 
performed for wild-type RRAS and the RASopathy causative 
mutant. The particle mesh Ewald method was used to evaluate 
electrostatic interactions (62), whereas a cut-off scheme was 
employed for Van der Waals interactions. Temperature and pres- 
sure were kept constant at 300 K and 1 bar by using the Berend- 
sen weak-coupling method (63), using separate temperature 
baths for protein-GDP complex and solvent, with a relaxation 
time of 0.1 ps for temperature and 1 ps for pressure. A time 
step of 2 fs was employed. Root mean square deviations were 
calculated according to standard definitions. UCSF Chimera 
(http://www.cgl.ucsf.edu/chimera/) was used for molecular 
graphics and structures superposition, by using the MatchMaker 
option. 

Biochemical and functional characterization of RRAS 
mutants 

The generation of constructs, and preparation and purification of 
proteins were as previously described (64). The intrinsic activ- 
ities of the RAS proteins, their modulation by GEFs and GAPs 
and their interaction with different effector proteins were deter- 
mined as described earlier (65,66). Dissociation of mantGDP 
from RAS proteins (0. 1 |jlm) was measured in the presence of 
20 |jlm GDP in 30 im Tris-HCl, pH 7.5, 10 mM KH 2 P0 4 / 
K 2 HP0 4 , 5 mM MgCl 2 , 3 mM dithioerythritol (DTE), at 25 °C, 
using a Perkin Elmer fluorimeter at 366 nm (excitation wave- 
length) and 450 nm (emission wavelength). Observed rate con- 
stants (ko bs ) of dissociation were obtained by single exponential 
fitting of the data. GEF-accelerated mantGDP dissociation from 
RAS proteins (0.1 |xm) was measured as mentioned earlier, in 
the presence of the catalytic domain of SOS1, Cdc25 (5 |jlm), 
using stopped-fiow instrument. 

The intrinsic GTPase reaction was performed by mixing 
70 |jlm nucleotide-free RAS proteins (HRAS, RRAS , 
RRAS V55M and RRAS G39dup ) with 50 pM GTP in 30 mM Tris- 
HCl, pH 7.5, 10 mM KH 2 P0 4 /K 2 HP0 4 , 10 mM MgCl 2 , 3 mM 
DTE, at 25 °C, using HPLC assay as previously described (67). 
Samples were taken at different time points and analysed by 
HPLC for their GDP and GTP contents to determine the relative 
GTP content [(GTP)/(GDP + GTP)]. Intrinsic GTP hydrolysis 
k 0 bs °f proteins were obtained by single exponential fitting of 
the data. For determination of GAP (neurofibromin, residues 
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l-333)-stimulated GTPase activity, GDP bound to HRAS and 
RRAS proteins was exchanged with excess mantGTP in the 
presence of alkaline phosphatase. Unbound nucleotides were 
removed by NAP5 column, and the RAS/mantGTP proteins 
were snap-frozen in liquid nitrogen (66). GAP-stimulated GTP 
hydrolysis of RAS proteins (0.2 jjlm) was measured in 30 mM 
Tris-HCl, pH 7.5, 10 mM MgCl 2 , 3 mM DTE at 25°C using a 
Hightech TgK Scientific stopped-fiow instrument. Reactions 
measured the decrease in fluorescence owing to hydrolysis of 
mantGTP. This decay was fit by a single exponential. 

Effector binding assays were performed in 30 mM Tris-HCl, 
pH 7.5, 100 mM NaCl, 5 mM MgCl 2 , 3 mM DTE at 25°C using a 
Fluoromax 4 fluorimeter in polarization mode. Increasing amounts 
of GST-tagged RAS-binding domains (RBD) of RAS effectors 
were titrated to 0.3 |jlm mantGppNHp-bound RAS proteins result- 
ing in an increase of polarization (64). The dissociation constants 
(IQ) were calculated by fitting the concentration-dependent 
binding curve using a quadratic ligand binding equation. 

For cell-based assays, COS-7 cells were transiently trans- 
fected with FLAG-tagged RRAS WT , RRAS V55M or RRAS G39dup 
by the DEAE-dextran method. For serum conditions, cells were 
incubated for 48 h in 10% FCS. In serum-starved conditions, 
serum was changed to basal medium midway between the trans- 
fection and harvesting. Transfected COS-7 cells were harvested 
and lysed in fishing buffer [50 mM Tris-HCl, pH 7.5, 2 mM 
MgCl 2 , 100 him NaCl, 1% IGEPAL CA-630, 10% glycerol, 
EDTA-free protease inhibitor cocktail (Roche, 1 tablet/50 ml 
buffer), 20 mM disodium B-glycerol phosphate and 1 mM 
Na 3 V04]. Cleared cell lysates were incubated with GSH-beads 
loaded with GST-RAF 1 -RBD. GTP-bound proteins and total re- 
combinant proteins were analysed by immunoblotting with anti- 
FLAG antibody. Antibodies against MEK1/2, ERK1/2, AKT, 
phospho-MEKl/2 (Ser217/221), phospho-ERKl/2 (Thr202/ 
Tyr204) and phospho-AKT (Thr308) were purchased from 
Cell Signaling Technology (68). 

Caenorhabditis elegans studies 

Culture and maintenance of animals were as previously 
described (69). The let-60(nl046) (let-60/RAS gain-of-function 
allele) and let-23(syl) (let-23/EGFR hypomorphic allele) 
strains were provided by the Caenorhabditis Genetics Center 
(University of Minnesota). The three-nucleotide insertion, 
c.81_82insGGC (ras-l G27dup ), corresponding to c.ll6_118dup 
in RRAS, was introduced in the wild-type cDNA (ras-l w ) 
(C. elegans ORF clone AAB03320, Thermo Scientific) by site- 
directed mutagenesis (QuikChange Site-Directed Mutagenesis 
Kit, Stratagene). ras-1 cDNAs were subcloned into the 
pPD49.83 heat shock-inducible vector (a gift of A. Fire, Stanford 
University School of Medicine). Germline transformation was 
performed as described (70). pJM371 plasmid [pelt-2:: 
NLSr.RFP] (a gift from J.D. McGhee, University of Calgary), 
which drives red fluorescent protein (RFP) expression in intes- 
tinal cell nuclei, was used as co-injection marker (30ng/|jil). 
Two different doses of constructs were injected (30 and 
100 ng/|xl). Animals from at least three independent transgenic 
lines for each construct and each dose of injection (i.e. six 
lines expressing ras-l WT and six lines expressing ras-l G27dup ) 
were heat-shocked in parallel and scored blindly at a Leica 
MZ10F dissecting microscope to check for the presence of 



protruding vulvae (Pvl phenotype) and multiple ectopic pseudo- 
vulvae (Muv phenotype), count the number of eggs retained in 
the uterus (Egl phenotype) and identify animals that had 
become bag-of-worms (Bag phenotype). Isogenic worms that 
had lost the transgene were cloned separately and used as con- 
trols. Following heat shock, all the transgenic lines expressing 
ras-1 WT or ras-l G27dup showed a variable degree of these pheno- 
types. Lines gbEx555a\hsp-16.41::ras-l WT ;pelt-2\-m,S::RFP^ 
and gbEx557a[hsp-16.41::ras-l G27dup ;pelt-2::NLS::RFP] were 
scored quantitatively in triplicate experiments at the compound 
microscope and used for further analyses. Genetic crosses were 
performed according to standard methods (69). The genotype of 
individual alleles was confirmed by direct sequencing of the ap- 
propriate genomic region. After each cross, isogenic worms that 
had lost the transgene were used as controls. 

To investigate VPCs induction and vulva morphogenesis, 
synchronized hermaphrodites carrying each transgene and the 
corresponding isogenic controls were heat-shocked in parallel 
at early L3 stage (33°C, 1 h, followed by 30°C, 1 h). Animals 
were scored at the compound microscope for vulval induction 
at late L3 and L4 stages, and for Pvl/Egl/Bag phenotypes at the 
adult stage. Microscopy observations were performed with a 
Nikon Eclipse 80i instrument equipped with Nomarski differen- 
tial interference contrast optics on live animals mounted on 2% 
agarose pads containing 10 mM sodium azide as anaesthetic. 

SUPPLEMENTARY MATERIAL 

Supplementary Material is available at HMG online. 
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